Chiral edge states and fractional charge separation in interacting bosons on a Kagome 
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We consider the extended hard-core Bose-Hubbard model on a Kagome lattice with boundary 
conditions on two edges. We find that the sharp edges lift the degeneracy and freeze the system into 
a striped order at 1 /3 and 2 /3 filling for zero hopping. At small hopping strengths, holes spontaneously 
appear and separate into fractional charges which move to the edges of the system. This leads to 
a novel edge liquid phase, which is characterized by fractional charges near the edges and a finite 
edge compressibility but no superfluid density. The compressibility is due to excitations on the 
edge which display a chrial symmetry breaking that is reminiscent of the quantum Hall effect and 
topological insulators. Large scale Monte Carlo simulations confirm the analytical considerations. 
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Frustrated systems are a rich playground in the search 
for new exotic phases and excitations, such as spin liquid 
phases, 1 '- Dirac strings in a spin-ice*^ and fractional 
charges in Kagome lattice antiferromagnets. 5 Recently, 
progress of indirect observations of fractional excitations 
has been made^ but their detection remains far from 
trivial since the controlled excitation and separation of 
fractional charges is difficult. We now show that by in- 
troducing sharp edges on two sides of a Kagome lattice 
with interacting bosons, fractional charges appear spon- 
taneously and are located close to the separate edges de- 
pending on their chirality. These chiral edge states are 
reminiscent of phenomena in quantum hall physics^ and 
topological insulators^ but on the Kagome lattice the 
topology and chirality is completely controlled by the de- 
sign of the edges. Moreover, the appearance of the chiral 
fractional charges with a deconfined interaction gives rise 
to a new compressible quantum edge phase. 

Exotic excitations in frustrated systems such as frac- 
tional charges and monopoles can often be understood 
from rather straight-forward geometrical arguments. 9 
They appear as local defects in the form of a rearrange- 
ment of real charges on the nontrivial backgroundj 4 ' 5 ' 10 ! 11 
When two (or more) fractional charges separate they may 
be connected by a string of a slightly disturbed quantum 
ground state, which normally acts confining. This poses 
a generic problem for the controlled excitation and ob- 
servation of the fractional charges: the confining string 
keeps them close together so that a separate observa- 
tion becomes impossible, analogously to the difficulty of 
detecting individual quarks. If the string can be tuned 
to become deconfining, then the system typically under- 
goes a quantum phase transitions to a complicated sea of 
closed strings with entirely different properties. 

We now consider a system of interacting bosons on a 
Kagome lattice where boundary conditions play an im- 
portant role due to the huge degeneracy, and deconfin- 
ing strings appear, which allow a controlled separation 
of the fractional charges. The model is defined by hard 



core bosons on a Kagome lattice with nearest neighbor 
repulsion V and chemical potential \i 
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where t is the hopping between boson creation and an- 
nihilation operators b\ and bj on nearest neighbor sites 
The Kagome lattice is special since it allows a 
macroscopic degeneracy even for commensurate fillings 
of !/3 and 2 /3. In particular, without hopping the model 
in Eq. ([T]) corresponds to the Ising model with an entropy 
per sitei^ of S = 0.108 at filling 1/3, since all configu- 
rations with one particle per triangle are ground states 
i.e. there is a "triangle rule" analogous to the "ice rule" 
in the spin ice^ In contrast, the Ising model on the tri- 
angular lattice, which is also frustrated, does not show 
a macroscopic degeneracy at commensurate filling. Fi- 
nite hopping t > lifts this degeneracy which leads to 
an ordered state with a finite structure factor. At even 
larger hopping a phase transition to a superfluid phase 
is observed, which is believed to be weakly first order or 
possibly second orderJ^ At negative hopping t = — v ji 
the model corresponds to the Heisenberg model which is 
at half filling one of the most promising candidate for a 
spin liquid^ 

Fractional charges on the Kagome lattice which are 
connected by strings have been extensively discussed^ for 
the fermionic model in Eq. (JTJ at filling 1 /3. Basically 
the fractional charges correspond to one empty triangle, 
which have two possible chiralities: up or down. Since 
one up and one down triangle together correspond to a 
single missing charge (see Fig. [TJj), it is clear that a sin- 
gle empty triangle has negative fractional bosonic charge 
-!/2. However, so far there is no proposal on how to di- 
rectly observe these exotic excitations. 

As we will show here, it is possible to make a con- 
trolled separation of fractional charges by introducing 
sharp edges in the model of Eq. (JTJ on two sides of the 
Kagome lattice, which also gives rise to a new boundary- 
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FIG. 1: (a) Striped solid phase induced by edges, (b) One 
hole corresponds to two fractional charges (black triangles), 
(c) Hopping of real particles (blue dots) separates the frac- 
tional charges and resonant hexagons appear along the con- 
necting string (red dots), (d) The spin (red arrows) represen- 
tation of the string (yellow path), (e) Resonant hexagons can 
be flipped by a third order hopping process. 

induced quantum edge phase. Interestingly, such edges 
as depicted in Fig. [1] have a macroscopic effect in the 
Kagome lattice, since the edges fix the bosons at the 
upper corners of the triangles. Therefore the order be- 
comes frozen for t = and there is a unique ground state 
of bosons in stripes parallel to the edges as shown in 
Fig. [TJ,. In other words, the macroscopic entropy is lost 
in this "striped solid phase" , even when the edges are 
very far apart. 

In a system with such edges the role of finite but small 
hopping t is now dramatically reversed: Instead of lifting 
the degeneracy, the hopping now facilitates fluctuations, 
which turn out to be nothing else but strings between 
fractional charges. As depicted in Fig. [JJ the separation 
of two fractional charges facilitates a string with configu- 
rations of triple occupied hexagons. At finite hopping the 
triple occupied hexagons are preferred since they lower 
the energy by — 12t 3 /V 2 via a resonant third order hop- 
ping process^ as shown in Fig. QJ. 

At this point two interesting observations can be made: 
First of all there is a new competition between a frozen 
ground state due to the edges without any resonant 
hexagons (striped solid phase) on the one hand and pos- 
sible fluctuating strings on the other hand. Secondly, 
closed strings cannot appear by a local rearrangement of 
charges without violating the triangle rule on the frozen 
configuration shown in Fig. QJi. Indeed it is straight- 
forward to convince oneself that fluctuating strings with 
triple occupied hexagons can only occur if either the tri- 
angle rule is violated or holes are introduced into the 
system. For the case /i < V we therefore expect that 
holes may appear spontaneously which separate into two 
fractional charges since the connecting strings allow ad- 
ditional resonant hexagons and therefore are deconfining. 
However, this phase is not the superfluid phase since the 



fractional charges are not free. Indeed we will see that 
the fractional charges only appear as edge excitations, so 
we will call this phase the "edge liquid phase" . 

In order to establish the new edge phase quantitatively, 
it is important to understand the process of fractional 
charge separation in detail as illustrated in Fig. [TJ In 
Fig. [Hi the frozen configuration in the striped order is 
shown, which is the unique ground state for t = 0. If 
a charge is removed (Fig. [TJd) two fractional charges can 
separate for t ^ 0, namely one up and one down trian- 
gle (Fig. [TJ;) which are connected by a string (red dots). 
Because of this string, the down triangle can never move 
underneath the up triangle. Since there are now reso- 
nant configurations (hexagons with three bosons) along 
the string, it is deconfining and the energy can be low- 
ered by pushing the fractional charges to the upper and 
lower edges respectively. When the energy gain from a 
fluctuating string exceeds the energy cost fj, of removing a 
charge, it can be expected that fractional charges appear 
spontaneously at the edges, leading to a new quantum 
edge phase which is compressible. 

The energy of a string can in turn be quantitatively 
estimated by realizing that in fact the resonances on each 
hexagon correspond to fluctuations of the string as shown 
in Fig. [TJi. Here the path of the string is mapped to up 
and down spins depending on if the string is to the right 
or left of a resonant hexagon. A flip can occur when 
the neighboring hexagon has opposite spin, leading to 
an xy— type model with effective exchange of — 12t 3 /V 2 
along the length of the string L s < L y 

10+3 L = 

H *y = — y^J2( S tSi +1 + h.c). (2) 

i—1 

For a given length L s the ground state energy of this 
model is given by 

, . 12t 3 ( 7T \ 24i 3 

E, y (L s ) = ^(l-co iw - T) ^-— Lai (3) 

which means that there can be a substantial energy gain 
by maximizing the length L s . However, this is only part 
of the story since the ends of the string can move freely up 
to the total systems size L y > L s , due to the hopping of 
the fractional charges. This hopping gives an additional 
kinetic energy of order t and independent of length, but 
for large systems sizes the leading effect of the xy— model 
in Eq. ^ is to maximize the length of the string L s , 
i.e. push the fractional charges apart towards the edges 
by an effective linear repulsive potential between the frac- 
tional charges. The interplay of maximizing L s and the 
kinetic energy of the fractional charges leads to a char- 
acteristic maximum of fractional charges near the edges, 
which resembles the solution of a particle in a linear po- 
tential (Airy function). At the top edge the fractional 
charges appear as down-triangles only and at the bot- 
tom edge they appear as up-triangles. These edge exci- 
tations therefore posses a definite chirality, which relates 
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FIG. 2: The phase diagram of the Bose-Hubbard model with 
cylindrical boundary conditions for L y = 18 and 24 from 
Monte Carlo simulations at temperature T — 0.005V. The 
phase transition line to the superfiuid phase was determined 
on lattices up to L y = 24 with periodic boundary conditions. 



the pseudo-spin of the fractional charge (up or down) to 
the direction given by the edge. 

We now turn to quantum Monte Carlo simulations of 
the hard-core boson model in Eq. (p} on the Kagome 
lattice in order to numerically analyze the existence of 
the edge excitations. We use the stochastic cluster se- 
ries expansio n 14 i 15 with parallel temperingi£~— on sys- 
tem sizes up to 1752 sites (L y = L x = 24) with cylin- 
drical boundary conditions. In Fig. [5] we show a small 
section of the phase diagram^ near the transition to 
the 1 /3 solid phase for cylindrical boundary conditions. 
As predicted, the striped solid order is destroyed by the 
spontaneous appearance of fractional charges as the hop- 
ping is increased. The value of the critical hopping is 
lower for larger L y , which is in agreement with the ar- 
gument that the energy of the string is proportional to 
L s . In fact, the effective model above can also be used to 
determine the critical hopping, which gives a compatible 
estimate at about 5% larger hopping. In the thermody- 
namic limit, the energy of strings will always be lower 
than the chemical potential and the edge liquid phase 
exists for all hopping strengths. 

In the middle panel of Fig. [3J a typical snapshot of 
a configuration during a Monte Carlo simulations for 
t = 0.102V in the edge liquid phase is shown. A string is 
clearly visible by the shaded resonant hexagons (green) 
and the shaded up and down triangles (red) which in- 
dicate fractional charges near the edges. The distribu- 
tion of up-triangle fractional charges is shown in the left 
panel of Fig. [3] with a characteristic maximum at the 
lower edge, which agrees well with the results from the 
effective model discussed in Eq. ([2]). In the striped solid 
phase (t = 0.09V) no fractional charges are found (ex- 
cept for virtual excitations) and in the superfiuid phase 
(t = 1.2V) there are fractional charges in the entire sam- 
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FIG. 3: A snapshot of a typical configuration during a Monte 
Carlo run at p/V = 0.8, t/V = 0.102 and /3 = 200 with 
L x = L y = 24 is shown in the middle panel. Red indi- 
cates a higher density of empty triangles or fractional charges 
and green represents a higher density of hexagons in a res- 
onant configuration. The density distribution of up-triangle 
fractional charges is shown in the left panel for the striped 
solid (t = 0.09V), edge liquid (t = 0.102V) and superfiuid 
(t = 0.12) phases. The right panel shows the corresponding 
densities of resonant hexagons. 



pie, but the maximum near the edge remains. The right 
panel of Fig. [3J shows the change in resonant hexagon 
density in the different phases. 

The relevant order parameters in the different phases 
are shown in Fig. |H The spontaneous appearance of 
additional charges at a critical hopping is clearly seen 
by the change of total hole density Ap (relative to the 
striped order phase). The structure factor S(Q)/N = 
i\J2k=i n ke iQ ' rk \ 2 )/N 2 at Q = (2tt, 0) is an indication 
of the striped order, which drops sharply when the frac- 
tional charges appear in the edge liquid phase, but re- 
mains finite. At still larger hopping the system enters 
the superfiuid phase, which has a spontaneously broken 
U(l) symmetry (off-diagonal order) and a finite super- 
fluid density (p a = (W 2 /2/3t) in terms of the winding 
number W in Monte Carlo simulations^). The com- 
pressibility per site k is maybe the most interesting order 
parameter, since it is zero in the striped order phase and 
then has a clear maximum close to the phase transition 
into the edge liquid phase. The bottom panel of Fig. @] 
illustrates how this maximum (i.e. the phase transition) 
moves to lower values of t as the system size increases. 
At the same time, the compressibility per site decreases 
proportional to 1/L y , which is a clear indication that the 
compressibility is due to an "edge compressibility" . In- 
deed according to the considerations above the bulk sites 
do not contribute to the compressibility since there are 
no fractional charges (see Fig. [3J . The phase transition is 
basically characterized by the filling of interacting strings 
which makes it second order. The system is most com- 
pressible at the phase transition, but additional strings 
are harder and harder to insert which leads to a charac- 
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FIG. 4: (a) The change of hole density Ap relative to the 
striped solid phase, the structure factor with Q = (2tt, 0) and 
the superfhud density p s at fi/V = 0.8 and ft — 200 with 
Lx = L y = 24 in different phases, (b) The compressibility at 
(j,/V = 0.8 for different sizes. 

tcristic maximum. The oscillations of the compressibility 
as a function of hopping t are a finite size effect, where 
the addition of an extra string leads to regular blockades 
before the next string can be added, analogously to the 
quantized charging energy of a small capacitor (Coulomb 
blockade). 

In summary, the results show that edges have a dra- 
matic effect for hard core bosons on a Kagome lattice 
and lead to a new edge liquid phase, which is charac- 
terized by a finite edge compressibility but no superfhud 
density. Quantum Monte Carlo simulations and analyti- 
cal arguments demonstrate that fractional excitations of 
up- and down-triangles can be separated in space and 
are localized close to the lower and upper edges respec- 



tively. The fractional charges are connected by quantum 
strings of resonant configurations, which act deconfining. 
Edge states with a finite compressibility and the locking 
of an internal quantum number with the direction of the 
edge (e.g. the chirality of the spin-momentum locking) 
are also famous characteristics of topological insulators 
and the quantum Hall effect. However, for the hard-core 
boson system studied here the non-trivial topology is not 
a property of the bulk, but is instead controlled by the 
design of the sharp edges. In fact, the study of different 
kind of edges or defects in the edges are promising future 
research topics. 

Recently, much experimental progress has been made 
in the realizations of related models with ultra-cold 
bosonic gases. In particular, it is now possible to create 
an artificial Kagome lattice for ultra-cold bosons. 20 More- 
over, a gas of Rydberg atoms was successfully loaded in 
a two-dimensional optical lattice, and the spatially or- 
dered structures induced by the short range repulsion 
interactions have been observed directly- 2 ^. In principle 
it should be possible to create sharp boundaries on such 
systems. 22 Together with recent advances in single-site 
addressability 2 - - — , the effects illustrated above therefore 
give a promising perspective of actually taking pictures 
of strings and localized fractional charges analogous to 
the snapshot in our simulations in Fig. [3] 
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